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An analytical expression for the disc power is derived based on an equivalent circuit in black 
hole (BH) magnetosphere with a mapping relation between the radial coordinate of the disc and 
that of unknown astrophysical load. It turns out that this disc power is comparable with two other 
disc powers derived in the Poynting flux and hydrodynamic regimes, respectively. In addition, the 
^ , relative importance of the disc power relative to the BZ power is discussed. It is shown that the BZ 

power is generally dominated by the disc power except some extreme cases. Furthermore, we show 
^-j' that the disc power derived in our model can be well fitted with the jet power of M87. 

& : 

PACS numbers: 97.60.Lf, 98.62.Mw, 98.62.Js 

r-» ! 

As is well known, jets exist in many astronomical cases, such as active galactic nuclei, quasars, and 
£SJ . young stellar objects. The association of jets with magnetized accretion discs or magnetized central 

J> ' objects (black holes or stars) is strongly supported by the recent observations of Hubble Space Telescope 

, (HST), Chandra, and VLBlJ 1,2 ' Different theoretical models have been proposed for acceleration and 

collimation of jets, which can be divided into two main regimes, the Poynting flux regime and the 
| hydromagnetic regime I 3-6 !. Both regimes are related to a poloidal magnetic field threading the disc, 

■ from which energy and angular momentum are extracted. In the Poynting flux regime, energy is 
| extracted in Poynting flux (i.e. purely electromagnetic energy), but in the form of magnetically driven 

f^S . material winds in the latter regime. 

""^5 ' Recently, Cao discussed the disc power in the hydrodynamic regime, and found that the disc power 

(henceforth Pmhd) is mainly determined by the strength and the configuration of the poloidal magnetic 
field. I 7 ' Lee discussed the disc power in the Poynting flux regime (henceforth P EM ), and found that 
the Poynting flux caused by a rotating magnetic field with Keplerian angular velocity can balance the 

■ energy and angular momentum conservation of a stationary accretion flowJ 8 ! In this Letter, we derive 
the disc power (henceforth P^m) by using an equivalent circuit in black hole (BH) magnetosphere, 
and compare P E l M with Pmhd and P E m-> respectively. Throughout this Letter the geometric units 
G = c = 1 are used. 
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In order to derive the expression of the disc power P EM we make the following assumptions: 



EM 

■ 1. The disc is thin, Keplerian, stable and perfectly conducting, which lies in the equatorial plane of 

the BH with the inner boundary being at the marginally stable orbit with the radius rj„ = r TOS J 9 l 

2. The configuration of the poloidal magnetic field of the BH-disc system is shown in Fig. 1. The 
poloidal magnetic field lines thread both the BH horizon and the surrounding disc, connecting the disc 
with the remote astrophysical load. The shape of each field line threading the disc at radius r has an 
analogy with that in the asymptotic jet model. ^ It is assumed that the asymptotic radius of each field 
line r' is smaller than its light radius Rl- The asymptotic radius r' is related to the disc radius r by 

r' = Xr (1) 

3. The astrophysical load is axisymmetric, being located evenly in a plane P with some height above 
the disc. The surface resistivity ctl of the load obeys the following relation, ctl = az&H = 47ra^, where 
an = 4-7T = 377 ohm is the surface resistivity of the BH horizon, and az = 1 is taken in calculations. 

4. The poloidal magnetic field on the disc varies as a power law as given in Refs. [11-13], 

B£«r n , (2) 
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FIG. 1: Configuration of poloidal magnetic field threading the BH and its surrounding disc. 
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FIG. 2: Equivalent circuit for magnetic extraction of energy from accretion disc around a rotating BH. 



where £ = r/r ms is the radial parameter of the disc defined in terms of the radius of inner edge of the 
disc. The field threading the BH is just a continuation of the field threading the disc, and the following 
relation is satisfied, ^ 



<;B P D (r ms ) . 



(3) 



where B P H and B P D (r ms ) are the magnetic fields at the horizon and the inner edge of the disc, respectively. 
Some authors I 14 - 15 ! argued that the magnetic field threading the BH should not be significantly stronger 
than the field threading the inner disc. In this Letter q — 3 is assumed in calculations. 

Based on MacDonald and Thome's work,! 16 ] we propose an equivalent circuit to calculate the disc 
power as shown in Fig. 2. The segments PS and QR represent two adjacent magnetic surfaces consisting 
of the field lines connecting the disc and the load, and the segments PQ and RS represent the disc 
surface and the load sandwiched by the two surfaces, respectively. The quantities AZa and Aed 
(the subscript "i" is omitted) are the resistance of the load and emf due to the rotation of the disc, 
respectively. The disc load is neglected in calculations due to its perfect conductivity. 

Based on assumptions and J3J) the load resistance AZa between the two adjacent magnetic surfaces 
can be written as 



AZ, 



dr' 
2tt? 



dr 



(4) 
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The following equations are used in deriving disc power, which are almost the same as given in deriving 
the BZ power and MC power in our previous work. I 12 ' 13 ! 

AP I E I M = {P) 2 AZ A , P=Ae D /AZ A , Ae D = (A* D /2Tr)n D , (5) 

where I p is the current in each loop, and A^u is the magnetic flux between the two adjacent magnetic 
surfaces, 

AV D = B p D 2n (vjp I Va) dr, (6) 

V / / 9— 7r/2 

The quantity Qd is the angular velocity of the disc at the place where the magnetic flux penetrates, 
and it reads 



The concerned Kerr metric coefficients are given as follows^ 17 ! 

w = (S/p) sin#, 

p 2 = r 2 + a 2 cos 2 9, 
A = r 2 + a 2 - 2Mr. 



Combining Eqs.(2) with (9), we have the poloidal magnetic field on the disc expressed by 

Bl = ,- 1 BU- n . (9) 

Incorporating Eqs. (5) — (9), we have 

AP*? M /P = f(a*,t,ri)(%, (10) 
where the function /(a*,£,n) is expressed by 

y 8 (l + £- 2 y- 4 a 2 + 2a 2 £- 3 Y~ 6 ) £~ 2n+3 
f(a*,£,n) = ^. (11) 

2a z , 2 (1 + d- 2 x™Ul - 2^ X ™1) (e /2 X 3 ms + a.) 

The quantity Pq is defined as 

Pa = B 2 H M 2 w B 2 (M/Mq) 2 x 6.59 x 10 28 erg • s _1 , (12) 

where B4 and M/Mq are the strength of B P H and the BH mass in the units of I0 4 gauss and one solar 
mass, respectively. Integrating Eq.(10) over the radial parameter £, we have the expression for the disc 
power as follows, 

PMm (a*, e, n)/P =Jf (a*, n) (13) 
Cao derived the disc power in the hydromagnetic regime as follows, ^ 



3 



L *=\J (B p D f[rQ D (r)} a f(a, 7j )rdr, (14) 



/(a, 7 ,) = ( 7 ,-l) 7 «( 7 |-l)- (tt+1)/2 , (15) 

where the function / (a, 7 j) depends on the self-similar index a > 1 and the Lorcntz factor of the jet 
7j. Since both rf^D (r) and / (a, 7 j) are always less than unity, the upper limit to the disc power Ld 
can be written as 

U < Pmhd = \! {B p D f r 2 n D (r) dr. (16) 
The expression of the poloidal magnetic field is expressed as follows, ^ 

B P D (r) oc (r/M )" 3/4 (M/M )~ 1/2 A~ 1 BE 1 ^ 2 , (17) 

where A, B and E are general relativistic correction factors given in Rcf. [9]. Combining Eq.(16) with 
Eq.(17), we obtain 

where K = (A^BE 1 ' 2 )^ I (A^BE 1 / 2 )^ 

Lee derived the disc power P I EM in the Poynting flux regime as follows, ^ 

P'em = 4{[-«o (rout)} - [-«o (r m )}} [B P D (r m ) r ln f . (19) 
Setting the inner and outer edge of the disc at ri„ = r rns and r out = ^r ms , we have 

P'em (a*,$/P = 4^ 2 X f„ s (£ f - EL) , (20) 

where — — uq (£r ms ) and Ej ns = — uq (r ms ) are the specific energies of the disc matter corresponding 
to r and r ms: respectively. ^ 

We are going to compare our disc power P^ M with Pmhd and P^m- To facilitate comparison between 
PJ/m and Pmhd, we take n = 3/4 in Eq.(13), which is the same as given in Eq.(17). By using Eqs.(13), 
(18) and (20) we have the curves of the three disc powers varying with the radial parameter £ for the 
given values of shown in Fig. 3. 

In order to compare the above disc powers more clearly, we define 771 = Pem/Pmhd and 772 = 
P^m I 'Pemi ano - have the contours of 771 and 772 with different values in £ — a* parameter spaces as 
shown in Fig. 4 (a) and (b), respectively. 

Inspecting Figs. 3 and 4, we have the following results: 

(1) If the BH spin is not extremely high, the three disc powers satisfy the relation, P^/m < Pem < 
Pmhd- However this order reverses to Pmhd < Pem < Pemi when the BH spin approaches unity, 

(2) The disc power P EM dominates over P^ M and Pmhd in two cases: (i) The region is very close 
to the inner edge of the disc, (ii) The BH spin is very high. 

Recently some authors argued that the BZ power was overestimated, and it is, in general, dominated 
by the electromagnetic power output of the inner region of the disc, provided that the poloidal magnetic 
field threading the BH does not differ significantly in strength from that threading the disc. I 14 ' 15 ' In 
this Letter we shall discuss the importance of the disc power relative to the BZ power. The BZ power 
has been derived also based an equivalent circuit in BH magnctosphere in our previous workJ 12 ' 13 ! 

ft'/fl^ nv'?'^ (2i » 

Jo 2 - (1 - q) sin 
4 




where q = y/l — a% is a function of the BH spin, and k = ftp /fin is the ratio of the angular velocity 
of the magnetic field lines to that of the horizon. The angular coordinate 9 varies from to 9l on the 
BH horizon as shown in Fig.l. In calculations k = 0.5 and 9l = 0.45-7T are assumed in Eq.(21). In order 
to compare the disc power with the BZ power clearly, we defined the ratio of the disc power to the BZ 
power as r) = PJ/m/Pbz- For the given values of power-law index n the contours of r] are shown in 
£ — a* parameter space in Fig. 5. 

From Fig. 5 we find that the disc power will, in general, dominate over the BZ power, except that the 
BH spin a* is extremely high. This result is consistent with those given in Rcfs.[14] and [15]. 

In our model, different magnetic flux surfaces anchor at different foot point radii r on the disc, 
which correspond to different light cylinder radii R^. It is reasonable that the asymptotic radius of 
the magnetic flux surface is not greater than its light cylinder radius, i.e., r' < Rl- The light cylinder 
radius satisfy the condition, 



Rl^d = 1. (22) 

Incorporating Eqs.(7) and (22), we obtain the location of the light cylinder radius of magnetic flux 
surface as follows, 
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FIG. 5: The contours of constant values of 77 in £ — a, parameter space for n — 3/4, n = 1, n = 4/3 in (a), (b) 
and (c), respectively. 



R L = M 



(23) 



Fendt and Memola defined jet expansion rate as the ratio of the asymptotic jet radius to the foot 
point jet radius (the 'disc radius'), and the structure of jet is well in accord with the observation of the 
M87 jet. I 10 ! The observations of the HST show that the radius of the M87 jet is about 100i? s , [181 where 
Rs = 2GM/c 2 is Schwarzschild radius. Taking the jet expansion rate A = 10 as given in Ref.[10],we 
infer that the foot point jet radius is no greater than lORs- 

By using Eq.(13), we have the disc power corresponding to different BH spin a* as listed in Table 
1, where £ is the upper limit of integration in Eq.(13). In calculations the mass of the central BH is 
3 x 10 9 A/ Q as given in Refs.[l] and [2], and the magnetic field on the BH horizon is assumed to be about 
3 x 10 2 gauss as estimated in Ref.[3]. It turns out that the values of the disc power P% M in Table 1 are 
consistent with the value of the M87 jet power. I 19 ! 

In this Letter the interaction between the BH and the surrounding disc is neglected for simplification. 
In fact the magnetic coupling! 12 ' 13, 2 °lof th e gfj -with its surrounding disc will result in a further decrease 
of the BZ power, and the domination of the disc power over the BZ power will be strengthened in this 
case. We shall discuss this situation in our future work. 



TABLE I: Values of 

Rem m accord with the M87 jet power for different values of the BH spin a*. 



a* 





0.3 


0.5 


0.75 


0.998 


Tms/Rs 


3 


2.489 


2.117 


1.579 


0.618 


£ = WR s /r ms 


3.333 


4.017 


4.725 


6.333 


16.169 


PUm^tq ■ s- 1 ) 


2.04 x 10 44 


1.91 x 10 44 


1.78 x 10 44 


1.57 x 10 44 


2.06 x 10 44 
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